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Abstract —We study the problem of efficiently broadcasting 
packets in multi-hop wireless networks. At each time slot the 
network controller activates a set of non-interfering links and 
forwards selected copies of packets on each activated link. A 
packet is considered jointly received only when all nodes in 
the network have obtained a copy of it. The maximum rate 
of jointly received packets is referred to as the broadcast 
capacity of the network. Existing policies achieve the broadcast 
capacity by balancing traffic over a set of spanning trees, 
which are difficult to maintain in a large and time-varying 
wireless network. We propose a new dynamic algorithm that 
achieves the broadcast capacity when the underlying network 
topology is a directed acyclic graph (DAG). This algorithm is 
decentralized, utilizes local queue-length information only and 
does not require the use of global topological structures such 
as spanning trees. The principal technical challenge inherent 
in the problem is the absence of work-conservation principle 
due to the duplication of packets, which renders traditional 
queuing modelling inapplicable. We overcome this difficulty by 
studying relative packet deficits and imposing in-order delivery 
constraints to every node in the network. Although in-order 
packet delivery, in general, leads to degraded throughput in 
graphs containing cycles, we show that it is throughput optimal 
in DAGs and can be exploited to simplify the design and 
analysis of optimal algorithms. Our characterization leads to a 
polynomial time algorithm for computing the broadcast capacity 
of any wireless DAG under the primary interference constraints. 
Additionally, we propose a multiclass extension of our algorithm 
which can be effectively used for broadcasting in any network 
with arbitrary topology. Simulation results show that the our 
algorithm has superior delay performance as compared to the 
tree-based approaches. 

I. Introduction and Related Work 

Broadcast refers to the fundamental network function¬ 
ality of delivering data from a source node to all other 
nodes. For efficient broadcasting, we need to use appropriate 
packet replication and forwarding to eliminate unnecessary 
packet retransmissions. This is especially important in power- 
constrained wireless systems which suffer from interference 
and collisions. Broadcast applications include mission-critical 
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military communications m , live video streaming m, and data 
dissemination in sensor networks 0. 

The design of efficient wireless broadcast algorithms faces 
several challenges. Wireless channels suffer from interfer¬ 
ence, and a broadcast policy needs to activate non-interfering 
links at every time slot. Wireless network topologies undergo 
frequent changes, so that packet forwarding decisions must 
be made in an adaptive fashion. Existing dynamic multicast 
algorithms that balance traffic over spanning trees i) may 
be used for broadcasting, since broadcast is a special case 
of multicast. These algorithms, however, are not suitable for 
wireless networks because enumerating all spanning trees 
is computationally prohibitive that needs to be performed 
repeatedly when the network topology changes with time. 

In this paper, we study the fundamental problem of through¬ 
put optimal broadcasting in wireless networks. We consider a 
time-slotted system. At every slot, a scheduler decides which 
non-interfering wireless links to activate and which set of 
packets to forward over the activated links, so that all nodes 
receive packets at a common rate. The broadcast capacity is 
the maximum common reception rate of distinct packets over 
all scheduling policies. To the best of our knowledge, there 
does not exist any capacity-achieving scheduling policy for 
wireless broadcast without the use of spanning trees The 
main contribution of this paper is to design provably optimal 
wireless broadcast algorithms that does not use spanning trees 
when the underlying topology is a DAG. 

We start out with considering a rich class of schedul¬ 
ing policies n that perform arbitrary link activations and 
packet forwarding. We define the broadcast capacity A* as 
the maximum common rate achievable over this policy class 
n. We next enforce two constraints that lead to a smaller 
set of policies. First, we consider the subclass of policies 
jjin-order ^ gnforce the in-order delivery of packets. 

Second, we focus on the subset of policies 11* C 
that allows the reception of a packet by a node only if 
all its incoming neighbours have received the packet. It is 
intuitively clear that the policies in the more structured class 
n* are easier to describe and analyze, but may yield degraded 
throughput performance. We show the surprising result that 
when the underlying network topology is a directed acyclic 
graph (DAG), there is a control policy tt* G H* that achieves 

’Note that we exclude network-coding operations throughout the paper. 
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the broadcast capacity. In contrast, we prove the existence of 
a network containing a cycle in which no control policy in the 
policy-space can achieve the broadcast capacity. 

To enable the design of the optimal broadcast policy, we es¬ 
tablish a queue-like dynamics for the system-state, represented 
by relative packet dehcits. This is non-trivial for the broadcast 
problem because explicit queueing structure is difficult to 
dehne in the network due to packet replications. We subse¬ 
quently show that, the problem of achieving the broadcast 
capacity reduces to hnding a scheduling policy stabilizing the 
system, which can be accomplished by stochastic Lyapunov 
drift analysis techniques IS), ||6|. 

In this paper, we make the following contributions: 

• We dehne the broadcast capacity of a wireless network 
and show that it is characterized by an edge-capacitated 
graph Q that arises from optimizing the time-averages 
of link activations. For integral-capacitated DAGs, the 
broadcast capacity is determined by the minimum in¬ 
degree of the graph Q, which is equal to the maximal 
number of edge-disjoint spanning trees. 

• We design a dynamic algorithm that utilizes local queue- 
length information to achieve the broadcast capacity of 
a wireless DAG network. This algorithm does not rely 
on spanning trees, has small computational complexity 
and is suitable for mobile networks with time-varying 
topology. This algorithm also yields a constructive proof 
of a version of Edmonds’ disjoint tree-packing theorem 
171 which is generalized to wireless activations but spe¬ 
cialized to DAG topology. 

• Based on our characterization of the broadcast capacity, 
we derive a polynomial-time algorithm to compute the 
broadcast capacity of any wireless DAG under primary 
interference constraints. 

• We propose a randomized multiclass extension of our 
algorithm, which can be effectively used to do broadcast 
on wireless networks with arbitrary underlying topology. 

• We demonstrate the superior delay performance of our 
DAG-policy, as compared to centralized tree-based algo¬ 
rithm ID, via numerical simulations. We also explore the 
efficiency/complexity trade-off of our proposed multiclass 
extension through extensive simulations. 

In the literature, a simple method for wireless broadcast is to 
use packet flooding na. The flooding approach, however, leads 
to redundant transmissions and collisions, known as broadcast 
storm M- In the wired domain, it has been shown that 
forwarding useful packets at random is optimal for broadcast 
this approach does not extend to the wireless setting due 
to interference and the need for scheduling im. Broadcast on 
wired networks can also be done using network coding na, 
ini- However, efficient link activation under network coding 
remains an open problem. 

The rest of the paper is organized as follows. Section |I^ 
introduces the wireless network model. In Section |In] we de¬ 
fine the broadcast capacity of a wireless network and provide 
a useful upper bound from a cut-set consideration. In Section 
IVl we propose a dynamic broadcast policy that achieves 


the broadcast capacity in a DAG. In section |V] we propose 
an efficient algorithm for computing the broadcast capacity 


of any wireless DAG under primary interference constraints. 
Our DAG-broadcast algorithm is extended to networks with 
arbitrary topology in section IVl] Illustrative simulation results 
are presented in Section |VIl| Finally, we conclude our paper 
in section Eul 

II. The Wireless Network Model 
We consider a time-slotted wireless network model repre¬ 
sented by the tuple (Q{V, E), c, , where V is the set of 
nodes, E is the set of directed links, c = (ce,e S E) is the 
capacity-vector of the links and S is the set of all feasible 
link-activation vectors, whose elements s = (se, e € E) € S 
are binary vectors such that the links e G E with Se = 1 can 
be activated simultaneously. The structure of the activation set 
S depends on the underlying interference model. For example, 
under the primary interference constraint (also known as node¬ 
exclusive interference constraint M), the set S consists of all 
binary vectors corresponding to matchings of the underlying 
graph g na, see Fig. [T] In the case of a wired network, S is 
the set of all binary vectors since there is no interference. 
In this paper we allow an arbitrary link-activation set S, 
which captures different wireless interference models. Let 
r S L be the source node at which stochastic broadcast traffic 
is generated (or arrives externally). The number of packets 
generated at the node r at slot t is given by the random variable 
A{t) G Z_|_, which is i.i.d. over slots with mean A. These 
packets are to be delivered efficiently to all other nodes in the 
network. 



(a) a wireless network 



(b) activation vector si 




(c) activation vector S 2 (d) activation vector S 3 


Fig. 1: A wireless network and its three feasible link activations 
under the primary interference constraint. 


III. Wireless Broadcast Capacity 
Intuitively, the network supports a broadcast rate A if there 
exists a scheduling policy under which all network nodes 
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can receive distinct packets at rate A. The broadcast capacity 
is the maximally supportable broadcast rate in the network. 
Formally, we consider a class 11 of scheduling policies where 
each policy tt S If consists of a sequence of actions {7rt}i>i 
executed at every slot t. Each action ttj comprises of two 
operations: (i) the scheduler activates a subset of links by 
choosing a feasible activation vector s{t) G 5; (ii) each node 
i forwards a subset of packets (possibly empty) to node j 
over an activated link (i,j) G l(s(f) = 1), subject to the link 
capacity constraint. The class If includes policies that use all 
past and future information, and may forward any subset of 
packets over a link. 

Let Rf{t) be the number of distinct packets received by 
node i G V from the beginning of time up to time t, under a 
policy TT G If. The time average liminf 7 ’_>,oo K^{T)/T is the 
rate of distinct packets received at node i. 

Definition 1. A policy tt is called a “broadcast policy of rate 
A” if all nodes receive distinct packets at rate A, i.e., 

minliminf = A, w. p. 1, (1) 

i&V T-^-oo T 

where A is the packet arrival rate at the source node r. 

Definition 2. The broadcast capacity A* of a wireless network 
is the supremum of all arrival rates A for which there exists 
a broadcast policy t: G IV of rate A. 

A. An upper bound on broadcast capacity A* 

We characterize the broadcast capacity A* of a wireless 
network by proving a useful upper bound. This upper bound is 
understood as a necessary cut-set bound of an associated edge- 
capacitated graph that reflects the time-averaged behaviour of 
the wireless network. We first give an intuitive explanation of 
the bound, assuming that the involved limits exist. Then in the 
proof of Theorem [T] we rigorously prove the bound by relaxing 
this assumption. 

Fix a policy tt G If. Let /3g be the fraction of time link 
e G E is activated under tt; that is, we dehne the vector 

1 ^ 

/3" = (/?:, eGE)= lim - ^ s^t), (2) 

t=l 

where s^{t) is the link-activation vector under policy tt in 
slot t. The average flow rate over a link e under the policy 
TT is upper bounded by the product of the link capacity and 
the fraction of time the link e is activated, i.e., Ce/3g. Hence, 
we can dehne an edge-capacitated graph Q'^ = (V, E, (Ce)) 
associated with policy tt, where each directed link e G E 
has capacity Cg = Ce/3g; see Fig. for an example of such 
an edge-capacitated graph. Next, we provide a bound on the 
broadcast capacity by maximizing the broadcast capacity on 
the ensemble of graphs Q'^ over all feasible vectors 

We dehne a proper cut U of the network graph as a 
proper subset of the node set V that contains the source node 
r. Dehne the link subset 

Eu = {{fj)GE\tGU,j^U}. (3) 

Since U C V, there exists a node n G V \ U. Consider the 
throughput of node n under policy tt. The max-how min-cut 



Fig. 2: The edge-capacitated graph for the wireless network 
with unit link capacities in Fig. [T] and under the time-average vector 
(3'^ = (1/2,1/4,1/4). The link weights are the capacities Ce/3/. The 
minimum proper cut in this graph has value 1/2 (when U = {r, a, c} 
or {r, b, c}). An upper bound on the broadcast capacity is obtained 
by maximizing this value over all vectors (3'^ G conv (5). 


theorem shows that the throughput of node n cannot exceed 
the total link capacity J^eeEu across the cut U. Since 

the achievable broadcast rate A^ of policy tt is an upper-bound 
on the throughput of all nodes, we have A’’^ < Eeerj. Cg/3-. 
This inequality holds for all proper cuts U and we have 

A’^ < min Ce/3g. (4) 

U: a proper cut ‘ ^ 

e^Eu 

Equation 0 holds for any policy tt G If. Thus, the broadcast 
capacity A* of the wireless network satishes 

A* = sup A’^ < sup min Cg j3f 

ttGU ti-GH U : a proper cut ^ 

e^Ejj 

< max min > Cg f3e, 

/3Gconv(5) U'. a proper cut 

eeEu 

where the last inequality holds because the vector as¬ 
sociated with any policy tt G H lies in the convex hull of 
the activation set S. Our hrst theorem formalizes the above 
intuitive characterization of the broadcast capacity A* of a 
wireless network. 

Theorem 1. The broadcast capacity A* of a wireless network 
Q (y, E, c) with activation set S is upper bounded as follows: 

X* < max ( min > Cg /3g I . (5) 

/3Gconv(5) \U: a. proper cut j 

e^Eu 

Proof of Theorem See Appendix ■ 

B. In-order packet delivery 

Studying the performance of any arbitrary broadcast policy 
TT G n is formidable because packets are replicated across 
the network and may be received out of order. To avoid 
unnecessary re-transmissions, the nodes must keep track of 
the identity of the received set of packets, which complicates 
the system state; because instead of the number of packets 
received, the system state is properly described here by the 
subset of packets received at each of the nodes. 

To simplify the system state, we focus on the subset of 
policies C H that enforce the following constraint: 
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Constraint 1 (In-order packet delivery). A network node is 
allowed to receive a packet p only if all packets {1,2,... ,p — 
1} have been received by that node. 


In-order packet delivery is useful in live media streaming 
applications El, where buffering out-of-order packets incurs 
increased delay that degrades video quality. In-order packet 
delivery greatly simplifies the network state space. Let Ri{f) 
be the number of distinct packets received by node i by time 
t. For policies in the set of received packets by time 

t at node i is {1,..., Therefore, the network state in 

slot t is given by the vector R{t) = [Ri[f),i G C). 

In section IV we will prove that there exists a throughput- 
optimal broadcast policy in the space when the un¬ 

derlying network topology is a DAG. Ironically, Lemma Q 
shows that there exists a network containing a cycle in which 
any broadcast policy in the space is not throughput 

optimal. Hence the space can not, in general, be ex¬ 

tended beyond DAGs while preserving throughput optimality. 


Lemma 1. Let X*in-order be the broadcast capacity of the policy 
subclass C n that enforces in-order packet delivery. 

There exists a network topology containing a directed cycle 

such that \* in-order < A*. 


Proof of Lemma See Appendix ■ 

We will return to the problem of broadcasting in networks 


with arbitrary topology in Section VI 


This is achieved by a broadcast policy that uses in-order 
packet delivery. In other words, we show that imposing the in- 
order packet delivery constraint does not reduce the broadcast 
capacity when the underlying topology is a DAG. 

From a computational point of view, the equality in Eqn. is 
attractive, because it implies that for computing the broadcast 
capacity of any wireless DAG, it is enough to consider only 
those cuts that separate a single (non-source) node from the 
source-side. Note that, there are only |V| — 1 of such cuts, in 
contrast with the total number of cuts, which is exponential in 
the size of the network. This fact will be exploited in section 
[V]to develop a strongly poly-time algorithm for computing the 
broadcast capacity of any DAG under the primary interference 
constraints. 


IV. DAG Broadcast Algorithm 
In this section we design an optimal broadcast policy for 
wireless DAGs. We start by imposing an additional constraint 
that leads to a new subclass of policies H* C jjin-order 
we will see, policies in H* can be described in terms of 
relative packet deficits which constitute a simple dynamics. We 
analyze the dynamics of the minimum relative packet deficit 
at each node j, where the minimization is over all incoming 
neighbours of j. This quantity plays the role of virtual queues 
in the system and we design a dynamic control policy that 
stabilizes them. The main result of this section is to show that 
this control policy achieves the broadcast capacity whenever 
the network topology is a DAG. 


C. Achieving the broadcast capacity in a DAG 

At this point we concentrate our attention to Directed 
Acyclic Graphs (DAGs). Graphs in this class are appealing 
for our analysis because they possess well-known topological 
ordering of the nodes ITSl . For DAGs, the upper bound Q 
on the broadcast capacity A* in Theorem [T] will be simplified 
further. For each receiver node v f r, consider the proper cut 
Uv that separates the network from node v: 

Uo = V\ {u}. (6) 

Using these cuts {Uv,v rf r}, we define another upper bound 
Adag on the broadcast capacity A* as: 

Adag = max min Ce Pe (7) 

/3Gconv(»S) {Uv,v^r} 

eGEu.^j 

> max min > Cg /5e > A*, 

/3Gconv(5) U: a proper cut 

sGEu 

where the first inequality uses the subset relation {Uy,v rf 
r} C {U: a proper cut} and the second inequality follows 
from Theorem In Section |IV[ we will propose a dynamic 
policy that belongs to the policy class Hin-ordei and achieves 
the broadcast rate Adag- Combining this result with Q, we 
establish that the broadcast capacity of a DAG is given by 

A* = Adag = max min Ce {3e, (8) 

/3Gconv(tS) {Uv,v^r} 

e^Eu^j 

= max min > Cg Pe- 

/3Sconv(S) U-. a proper cut 

eCiEu 


A. System-state by means of packet deficits 


We showed in Section III-B that, constrained to the policy- 
space the system-state is completely represented by 

the vector R{f). To simplify the system dynamics further, we 
restrict further as follows. 

We say that node i is an in-neighbor of node j iff there 
exists a directed link {i,j) G E in the underlying graph Q. 


Constraint 2. A packet p is eligible for transmission to node 
j at a slot t only if all the in-neighbours of j have received 
packet p in some previous slot. 

We denote this new policy-class by H* C We will 

soon show that it contains an optimal policy. Fig. shows the 
relationship among different policy classe^ 

Following properties of the system-states Rit) under a policy 
TT € n* will be useful. 


Lemma 2. For j f r, let In{j) denote the set of in-neighbors 
of a node j in the network. Under any policy tt G H*, we 
have: 

(1) Rj{t) < minig/„(j) Rfit) 

(2) The indices of packets p that are eligible to be transmit¬ 
ted to the node j at slot t is given by 

{p \ Rj(t) -\- 1 < p < min Ri(t )}. 

iGIn(j) 


^We note that, if the network contains a directed cycle, then a deadlock 
might occur under a policy in H* and may yield zero broadcast throughput. 
However, this problem does not arise when the underlying topology is a DAG. 
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IT: all policies that perform 
link activations and routing 

j^in-ordt.r, enforce 

in-order packet delivery 

n*: policies that allow reception 
only if all in-neighbors have 
received the specific packet 


Fig. 3: Containment relationships among different policy classes. 

We define the packet deficit over a directed link {i,j) G E 
by Qij(t) = Ri{t) — Rj{t). Under a policy in If*, Qifit) 
is always non-negative because, by part (1) of Lemma we 
have 

Qij{t) = Ri{t) - Rj{t) > min Rk{t) - Rj{t) > 0. 

fcein(i) 

The quantity Qij (t) denotes the number of packets received by 
node i but not by node j, upto time t. Intuitively, if all packet 
deficits Qij{t) are bounded asymptotically, the total number 
of packets received by any node is not lagging far from the 
total number of packets generated at the source; hence, the 
broadcast throughput will be equal to the packet generation 
rate. 

To analyze the system dynamics under a policy in H*, it is 
useful to define the minimum packet deficit at node j r hy 

Xj{t) = min Qij{t). (9) 

iGIn(y) 

From part (2) of Lemma Xj (t) is the maximum number of 
packets that node j is allowed to receive from its in-neighbors 
at slot t. As an example. Fig. shows that the packet deficits at 
node j, relative to the upstream nodes a, b, and c, are Qaj (f) = 
8, Qbj (t) = 5, and Q^j (t) = 4, respectively. Thus Xj (t) = 4 
and node j is only allowed to receive four packets in slot t 
due to Constraint We can rewrite Xj (t) as 

Xj{t) = Qi»j{t), where *(= arg min Qij{t), (10) 

* zGIn(j) 

and the node ij" is the in-neighbor of node j from which node 
j has the smallest packet deficit in slot f; ties are broken 
arbitrarily in deciding Our optimal broadcast policy will 
be described in terms of the minimum packet deficits {Xj{t)}. 



Fig. 4; Under a policy tt G 11*, the set of packets available 
for transmission to node j in slot t is {11,12,13,14}, which 
are available at all in-neighbors of node j. The in-neighbor of 
j inducing the smallest packet deficit is = c, and Xj{t) = 
mm{Qaj{t),Qbj{t),Qcj{t)} = 4. 


^We note that the minimizer i* is a function of the node j and the time 
slot t; we slightly abuse the notation by neglecting j to avoid clutter. 


B. The dynamics of the system variable Xj (t) 

We now analyze the dynamics of the system variables 

X,{t) = Q.^^,(t) = ReSt)-RAt) (11) 

under a policy tt G 11*. Define the service rate vector /r(f) = 

iTij(f)){i,j)eE by 

Cij if {i,j) G E and the link {i,j) is activated, 

0 otherwise. 

Equivalently, we may write p,ij{t) = CijSijit), and the number 
of packets forwarded over a link is constrained by the choice 
of the link-activation vector s{t). At node j, the increase in the 
value of Rj{t) depends on the identity of the received packets; 
in particular, node j must receive distinct packets. Next, we 
clarify which packets are to be received by node j at time t. 

The number of available packets for reception at node 
j is ^i'o.{Xj{t),J2kev because: (i) Xj(t) is the 

maximum number of packets node j can receive from its in¬ 
neighbours subject to the Constraint El (11) Is 

the total incoming transmission rate at node j under a given 
link-activation decision. To correctly derive the dynamics of 
Rj{t), we consider the following efficiency requirement on 
policies in 11*: 

Constraint 3 (Efficient forwarding). Given a service rate 
vector p(t), node j pulls from the activated incoming links 
the following subset of packets (denoted by their indices) 

\ Rjff) + 1 < p < Rj{f) + min{Ai^(f), E 

kev 

( 12 ) 

The specific subset of packets that are pulled over each 
incoming link are disjoint but otherwise arbitrary^ 

Constraint requires that scheduling policies must avoid 
forwarding the same packet to a node over two different 
incoming links. Under certain interference models such as 
the primary interference model, at most one incoming link 
is activated at a node in a slot and Constraint [3 is redundant. 

In Eqn. the packet deficit Qi^jf) increases with 

Ri*{t) and decreases with Rj{t), where Ri*{f) and Rj{f) 
are both non-decreasing. Hence, we can upper-bound the 
increment of Qijj{f) by the total capacity ktmi'it) 

of the activated incoming links at node Also, we can 
express the decrement of Qi*j{t) by the exact number of 
distinct packets received by node j from its in-neighbours, 
and it is given by min{Arj(f), 'f2k^v Tkj{t)} by Constraint 
Consequently, the one-slot evolution of the variable Qijj{f) is 
given b>|!] 

1 ) — ~ ^ ^ Tkj(t)) + 'y ^ 

kGV m£V 

= (^Xj[t) — 'y ^ Bkj(f)) + y ^ Trail if) J 

fcGV rnGV 

(13) 

^Due to Constraints and the packets in iig have been received by all 
in-neighbors of node j. 

^We emphasize that the node is defined in depends on the particular 
node j and time t, and may be different from the node 
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where (a;)'*' = inax{a;, 0) and we recall that Xj{t) — 

It follows that Xj{t) evolves over slot t according to 

(a) ib) 

Xj{t + 1) = min Qij{t + 1) < + 1) 

iGIn(j) 

+ E ( 14 ) 

kev mev 

where the equality (a) follows the definition of Xj{t), equality 
(b) follows because node i’l G In(j) and equality (c) follows 
from Eqn. ([T3j. In Eqn. •El, if il = r, we abuse the notation 
to define = ^{t) for the source node r, where 

A{t) is the number of exogenous packet generated at slot t. 

C. The optimal broadcast policy 

Our broadcast policy is designed to keep the minimum 
deficit process X{t) stable. Eor this, we regard the variables 
Xj{t) as virtual queues that follow the dynamics ([T4]i. By 
performing drift analysis on the virtual queues Xj(t), we 
propose the following max-weight-type broadcast policy tt*, 
described in Algorithmic We have tt* G II* and it enforces 
the constraints [C |C We will show that this policy 

achieves the broadcast capacity A* of a wireless network over 
the general policy class If when the underlying topology is a 
DAG. 


Algorithm 1 Optimal Broadcast Policy tt* for a Wireless 
DAG:_ 

At each slot t, the network-controller observes the state- 
variables G V} and executes the following ac¬ 

tions 

1: Eor each link {i,j) G E, compute the deficit Qij{t) = 
Ri{t) — Rj{t) and the set of nodes Kj{t) C out(j) for 
which node j is their deficit minimizer, given as follows 

Kj{t) ^ {k gV \ j = aig min Qmkit)}. (15) 

mGln(fc) 

The ties are broken arbitrarily (e.g., in favor of the highest 
indexed node) in finding the argmin( ) in Eqn.lfTSll. 

2: Compute Xj{t) = minigin(j) Qijit) for j ^ r and assign 
to link {i,j) the weight 

W,,{t) ^ {X,{t) - E (16) 

keK,{t) 

3: In slot t, choose the link-activation vector s{t) = 
{Se{t),e G E) such that 

s(f) G argmax CeSeWe{t). (17) 

s^S 

e£E 


Step 1 



Q.a(t) =7 

KAt) 

= {a} 

Q.bM = 7 

KAt) 

= {b, c}* 

Qr-ait) = 8 

Kbit) 

= {0} 

Qab{t) = 0 

KAt) 

= {0} 

Qaa(t) = 1 



Qba{t) = 1 





= (X,(t) - Xb(t) - X,{t))+ = 6 
X5(t) = 0 Wrb(t) = {Xi(t))+ =Q 

WrAt) = (2fc(t))+ = 1 

= (W(t))+ = 0 
tVac(t) = (Xe(t))+ = 1 
Wba<,t) = (X,(t))+ = 1 


JXe(t) = 1 


Step 3 



Eig. 5: Running the optimal broadcast policy tt* in slot t in a wireless 
network with unit-capacity links and under the primary interference 
constraint. Step 1: computing the deficits Qijit) and a tie is 

broken in choosing node a as the in-neighbor deficit minimizer for 
node c, hence c G Ka{t); node b is also a deficit minimizer for node 
c. Step 2: computing Xj{t) for J 7^ r and Wij{t). Step 3: finding the 
link activation vector that is a maximizer in GD and forwarding the 
next in-order packets over the activated links. Step 4: a new packet 
arrives at the source node r and the values of {Rr{t + 1), Ra{t + 
1), Rb{t + 1), Rc{t + 1)} are updated. 


4: Every node j ^ r uses activated incoming 
links to pull packets {Rjit) + l,...,Rj{t) + 
Xj{t)}} from its in-neighbors 
according to Constraint]^ 

5: The vector {Rj{t),j G V) is updated as follows: 


We illustrate the Algorithm in an example in Eig. 
The next theorem demonstrates the optimality of the broadcast 
policy TT*. 


Rj{t+1) ^ 


Rj{t) + A(t), 

Rj{t) + minlXii CijS^j{t),Xj{t)}, 


j = r, 
j + r, 
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Theorem 2. If the underlying network graph Q is a DAG, then 
for any exogenous packet arrival rate A < \dac> the broadcast 
policy TT* yields 

■ V ■ rRf{T) ^ 

minlimini- = A, w.p. 1, 

i^V T-i-oo T 

where Xdac the upper bound on the broadcast capacity A* 
in the general policy class 11, as shown in 0- Consequently, 
the broadcast policy tt* achieves the broadcast capacity A* 
for any Directed Acyclic Graphs. 

Proof of Theorem ^ See Appendix ■ 


D. Number of disjoint spanning trees in a DAG 

Theorem [^provides an interesting combinatorial result that 
relates the number of disjoint spanning trees in a DAG to the 
in-degrees of its nodes. 

Lemma 3. Consider a directed acyclic graph G = {V, E) that 
is rooted at a node r, has unit-capacity links, and possibly 
contains parallel edges. The maximum number k* of edge- 
disjoint spanning trees in G is given by 

k* = min din(v), 

where di„{v) denotes the in-degree of the node v. 

Proof of Lemma See Appendix ■ 


V. Efficient Algorithm for Computing the 
Broadcast Capacity of a DAG 

In this section we exploit Eqn. 0 and develop an LP to 
compute the broadcast capacity of any wireless DAG network 
under the primary interference constraints. Although this LP 
has exponentially many constraints, using a well-known sepa¬ 
ration oracle, it can be solved in strongly polynomial time via 
the ellipsoid algorithm ifThll . 

Under the primary interference constraint, the set of feasible 
activations of the graphs are matchings 051. Eor a subset of 
edges E' C E, let ^ {Oj 1}'^' where (e) = 1 if 
e G E' and is zero otherwise. Let us define 


'Pmatchmg(0) = convexhull({x'“ |M is a matching in G}) (18) 
We have the following classical result by Edmonds 03 - 


Theorem 3. The set VmatchingiG) is characterized by the set 
of all 13 G such that : 


I3e > 

E s 

e^E[U] 


0 Ve G E; (19) 

1 Vu G U 

~ S UcV,\U\odd 


Here E\U] is the set of edge (ignoring their directions) with 
both end points in U, Si„(u) (Sout(u)) denotes the set of all 
incoming (outgoing) edges to (from) the vertex u G V. 

Hence following Eqn. 0. the broadcast capacity of a DAG 
can be obtained by the following LP : 


max A 

(20) 

Subject to. 

A < ^ CePe Vu G U \ {r} 

(21) 


/3 € 7^matching {G) 

(22) 


Prom the equivalence of optimization and separation (via 
the ellipsoid method), it follows that the above LP is poly¬ 
time solvable if there exists an efficient separator oracle for the 
constraints Since there are only linearly many con¬ 

straints (|U| — 1, to be precise) in ( [2l] i, the above requirement 
reduces to an efficient separator for the matching polytope 
( |22l ). We refer to a classic result from the combinatorial- 
optimization literature which shows the existence of such 
efficient separator for the matching polytope 

Theorem 4. n7\t There exists a strongly poly-time algorithm, 
that given Q = (ViE) and (3 : E ^ determines if (3 
satisfies or outputs an inequality from that is violated 
by (3. 

This directly leads to the following theorem. 

Theorem 5. There exists a strongly poly-time algorithm to 
compute the broadcast capacity of any wireless DAG under 
the primary interference constraints. 

The following corollary implies that, although there are 
exponentially many matchings in a DAG, to achieve the broad¬ 
cast capacity, randomly activating (with appropriate probabil¬ 
ities) only |£'| -P 1 matchings suffice. 

Corollary 1. The optimal broadcast capacity A* in a wireless 
DAG, under the primary interference constraints, can be 
achieved by randomly activating (with positive probability) at 
most jE^I -P 1 matchings. 

Proof: Let {X*,(3*) be an optimal solution of the 
LP ( [20| ). Hence we have (3* G 7^matchmg((/) = 
convexhull({x^^|M is a matching in G}). Since the poly¬ 
tope 7^inatching((/) IS a subset of by Caratheodory’s 

theorem nsi, the vector /3* can be expressed as a convex com¬ 
bination of at most |E^|-Pl vertices of the poly tope 7^matching(^/)^ 
which are matchings of the graph Q. This concludes the proof. 


VI. Broadcasting on Networks with Arbitrary 
Topology 

In this section we extend the broadcast policy for a DAG 
to networks containing cycles. Prom the negative result of 
Lemma [2 we know that any policy ensuring in-order packet 
delivery at every node cannot, in general, achieve the broadcast 
capacity of a network containing cycles. To get around this 
difficulty, we introduce the concept of broadcasting using 
multiple classes 1C of packets. The idea is as follows; each 
class k G 1C has a one-to-one correspondence with a specific 





permutation of the nodes; for an edge (a, 6) G E if the 
node a appears prior to the node b in the permutation 
(we denote this condition by a -<k h), then the edge (a, h) is 
included in the class k, otherwise the edge (a, b) ignored by 
the class k. The set of all edges included in the class k is 
denoted by C E. It is clear that each class k corresponds 
to a unique embedded DAG topology Q^iV^E^), which is a 
subgraph of the underlying graph Q{V,E). 

A new incoming packet arriving at the source node is admitted 
to some class k G 1C, according to some policy. All packets in a 
given class k G K, ctre broadcasted while maintaining in-order 
delivery property within the class k, however packets from 
different classes do not need to respect this constraint. Hence 
the resulting policy does not belong to the class H* in but 
rather to the general class H. This new policy keeps the best 
of both worlds: (a) its description-complexity is Q{kN), where 
for each class we essentially have the same representations as 
in the in-order delivery constrained policies and (b) by relaxing 
the inter-class in-order delivery constraint it has the potential 
to achieve the full broadcast capacity of the underlying graph. 
Hence the broadcast problem reduces to construction of mul¬ 
tiple classes (which are permutations of the vertices V) out 
of the given directed graph such that it covers the graph effi¬ 
ciently, from a broadcast-capacity point of view. In Algorithm¬ 
ic we choose the permutations uniformly at random with the 
condition that the source always appears at the hrst position 
of the permutation. 

Theorem 6. The multiclass broadcast Algorithm^ with K 
classes supports a broadcast rate of 

K 

= max miny^ (26) 

where we use the convention that (3^^ = 0 if (i,j) ^ E^^\ 

The right hand side of Eqn. ( |26] l can be understood as 
follows. Consider a feasible stationary activation policy ttstat 
which activates class I on the edge {i,j) f3ij fraction of time. 
Since, by construction, each of the class follows a DAG, 
lemma © implies that the resulting averaged graph has a 
broadcast capacity of X^ = miuj for the class k. 

Thus the total broadcast rate achievable by this scheme is 
simply X^ = J2k=i = Efe minj E* Given the K 

classes, following the same line of argument as in ( |20l i, we can 
develop a similar LP to compute the broadcast capacity ( |26l ) 
of all these fc-classes taken together in strongly poly-time. 

The proof of Theorem follows along the exact same 
line of argument as in Theorem (|^, where we now work with 
the following Lyapunov function L{Q{t)), which takes into 
account all k classes: 

K 

= (27) 

k=l j^r 

We then compare the multiclass broadcast algorithm with 
the stationary activation policy ttsxat above to show that the 
Multiclass broadcast algorithm is stable under all arrival rates 
below A. The details are omitted for brevity. 


Algorithm 2 Multiclass Broadcast Algorithm for General 
Topology 

Require: Graph Q{V,E), total number of classes K 

1: Generate K permutations of the nodes V uni¬ 

formly at random (with the source {r} at the hrst po¬ 
sition) and obtain the induced DAGs G^iy,E^), where 
e = (a, b) G E'^ iff a -<k b. 

2: For each permutation ^k, maintain a class k and the 
packet-counter variables at every node i = 

1,2,...,\V\. 

3: Each class observes intra-class packet forwarding con¬ 
straints (1), (2) and (3) described in sections |n^ and |iy] 

4: Dehne the state variables {Q'^{t), X^{t)} and compute 
the weights {W^(f)}, for each class k = 
exactly as in Eqn. ( [Thl l, where each class k considers the 
edges E^ only for Eqns. (HI) and ( fTh] ). 

5: An incoming packet to source r at time t joins the class 
I corresponding to 

E ^j(t) (23) 

jeK^t) 

6 : The overall weight for an edge e (taken across all the 
classes) is computed as 

We(f) = max W^(t) (24) 

k:eGE>‘ 

1: Activate the edges corresponding to the max-weight acti¬ 
vation, i.e., 

s{f) G argmaxy^ CeSeWglf). (25) 

sGS ^^ 
e^E 

8 : For each activated edge e G s(t), forward packets corre¬ 
sponding to a class achieving the maximum in Eqn. ([24li. 


Since the broadcast-rate X^ achievable by a collection of K 
embedded DAGs in a graph Q is always upper-bounded by 
the actual broadcast capacity A* of Q, we have the following 
interesting combinatorial result from Theorem (j^ 

Corollary 2. Consider a wired network, represented by the 
graph Q (V, E). For a given integer K > 1, consider K classes 
as in Theorem with being their corresponding 

edge-sets. Then, for any set of non-negative vectors {/3^}Ei 
with the following lower-bound for the 

broadcast capacity X* holds: 

K 

^*^E“i?EGj4 (28) 

k—1 i 

where we use the convention that /3^j = 0 if {i,j) ^ E^. 

The above corollary may be contrasted with Eqn. 0, which 
provides an upper bound to the broadcast capacity A*. 

VII. Simulation Results 

We present a number of simulation results concerning the 
delay performance of the optimal broadcast policy tt* in wire¬ 
less DAG networks with different topologies. For simplicity. 
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we assume primary interference constraints throughout this 
section. Delay for a packet is defined as the number of slots 
required for it to reach all nodes in the network, after its arrival 
to the source r. 


Diamond topology 

We first consider a 4-node diamond topology as shown 
Fig. Link capacities are shown along with the links. The 
broadcast capacity A* of the network is upper bounded by 
the maximum throughput of node c, which is 1 because at 
most one of its incoming links can be activated at any time. 
To show that the broadcast capacity is indeed A* = 1, we 
consider the three spanning trees {7i,72,73} rooted at the 
source node r. By finding the optimal time-sharing of all 
feasible link activations over a subset of spanning trees using 
linear programming, we can show that the maximum broadcast 
throughput using only the spanning tree 71 is 3/4. The 
maximum broadcast throughput over the two trees {71,71} is 
6/7, and that over all three trees {71,72,73} is 1. Thus, the 
upper bound is achieved and the broadcast capacity is A* = 1. 

We compare our broadcast policy tt* with the tree-based 
policy TTtree in a. While the policy TTtree is originally proposed 
to transmit multicast traffic in a wired network by balancing 
traffic over multiple trees, we slightly modify the policy TTtree 
for broadcasting packets over spanning trees in the wireless 
setting; link activations are chosen according to the max- 
weight procedure. See Fig. 5(a) for a comparison of the 
average delay performance under the policy tt* and the tree- 
based policy TTtree ovcf different subset of trees. The simulation 
duration is 10® slots. We observe that the policy tt* achieves 
the broadcast capacity A* = 1 and is throughput optimal. 


Mesh topology 

The broadcast policy tt* does not rely on the limited tree 
structures and therefore has the potential to exploit all degrees 
of freedom in packet forwarding in the network; such freedom 
may lead to better delay performance as compared to the tree- 
based policy. To observe this effect, we consider the 10-node 
DAG network subject to the primary interference constraint in 
Fig. For every pair of node {i,/}, 1 < * < J < 10, the 
network has a directed link from i to j with capacity (10 — i). 
By induction, we can calculate the number of spanning trees 
rooted at the source node 1 to be 9! « 3.6 x 10®. We choose 
five arbitrary spanning trees {%, 1 < * < 5}, over which the 
tree-based algorithm TTtree is simulated. Table demonstrates 
the superior delay performance of the broadcast policy tt*, 
as compared to that of the tree-based algorithm TTtree over 
different subsets of the spanning trees. It also shows that a 
tree-based algorithm that does not use enough trees would 
result in degraded throughput. 


Multiclass Simulation for Arbitrary Topology 

We randomly generate an ensemble of 500 wired networks 
(not necessarily DAGs), each consisting of A^ = 10 nodes and 
unit capacity links. By solving the LP corresponding to Eqn. 
(|26ll, we compute the fraction of the total broadcast capacity 



(c) Tree T 2 (d) Tree Ta 

Fig. 6: A wireless DAG network and its three embedded spanning 
trees. 



Fig. 7: Average delay performance of the optimal broadcast policy 
TT* and the tree-based policy Tr,ree that balances traffic over different 
subsets of spanning trees. 


achievable using K randomly chosen classes by the Multiclass 
Algorithm]^ of section VI The result is presented in Figure]^ 
It follows that a sizeable fraction of the optimal capacity may 
be achieved by using a moderate number of classes. However 
the number of required classes for achieving a certain fraction 
of the capacity increases as the broadcast capacity increases. 
This is because of the fact that increased broadcast capacity 
would warrant an increased number of DAGs to cover the 
graph efficiently. 


VIII. Conclusion 

We characterize the broadcast capacity of a wireless network 
under general interference constraints. When the underlying 
network topology is a DAG, we propose a dynamic algorithm 
that achieves the wireless broadcast capacity. Our novel design, 
based on packet deficits and the in-order packet delivery 
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A 

tree-based policy TTtree over the spanning trees: 

Ti Ti^Ts Ti^Ta Ti^Ta 

broadcast 
policy TT* 

0.5 

12.90 

12.72 

13.53 

16.14 

16.2 

11.90 

0.9 

1.3 X 10^ 

176.65 

106.67 

34.33 

28.31 

12.93 

1.9 

3.31 X 10^ 

1.12 X 10“^ 

4.92 X lO^ 

171.56 

95.76 

14.67 

2.3 

3.63 X 10^ 

1.89 X 10^ 

1.40 X 10^ 

1.76 X lO^ 

143.68 

17.35 

2.7 

3.87 X 10^ 

2.45 X 10^ 

2.03 X 10^ 

1.1 X 10^ 

1551.3 

20.08 

3.1 

4.03 X 10^ 

2.86 X 10^ 

2.51 X 10^ 

1.78 X 10^ 

9788.1 

50.39 


TABLE I: Average delay performance of the tree-based policy TTtree over different subsets of spanning trees and the optimal broadcast policy 

* 

TT . 







Fig. 8: The 10-node wireless DAG network and a subset of spanning 
trees. 


constraint, is promising for application to other systems with 
packet replicas, such as multicasting and caching systems. 
Future work involves the study of arbitrary networks, where 
optimal policies must be sought in the class If \ 
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Appendix 

A. Proof of Theorem 

Fix an e > 0. Consider a policy tt G 11 that achieves a broad¬ 
cast rate of at least X* — e defined in ([T]); this policy tt exists 
by the definition of the broadcast capacity A* in Definition]^ 
Consider any proper cut U of the network Q. By definition, 
there exists a node i ^ U. Let s'^{t) = {s'^(t),e G E) be 
the link-activation vector chosen by policy tt in slot t. The 
maximum number of packets that can be transmitted across 
the cut U in slot t is at most which is the 

total capacity of all activated links across U, and the link 
subset Eu is given in The number of distinct packets 
received by a node i by time T is upper bounded by the total 
available capacity across the cut U up to time T, subject to 
link-activation decisions of policy tt. That is, we have 

T T 

RUT) < E E • E ^"(9, (29) 

t—1 cGEu t—1 

where we define the vector u = (u^, e G E), = Cel[eG£;u]’ 

and a ■ b is the inner product of two vectors]^ Dividing both 
sides by T yields 


R^T) 

T 



It follows that 


(a) , , , 

A* — e < min lim inf 

T—¥oo 


R]{T) 

T 


< lim inf 
T —^oc 


RJjT) 

T 


< lim inf u ■ 
T—foo 



(30) 


where (a) follows that tt is a broadcast policy of rate at least 
A* — e. Since the above holds for any proper-cut u G U, we 
have 


A* — e < minliminf It • f —'Y (31) 

uGU T—^oo \ E ^ / 

\ t=l ' 

Now consider the following lemma. 

Lemma 4. For any policy tt G 11, there exists a vector (3'^ G 
conv[S) such that 


min lim inf M- f — | =mmu-l3^ w.p.l 

uClU T—foo \R ^^ / uClU 

t=l ■' 


Proof: Consider the sequence = L in¬ 

dexed by T > 1. Since s^{t) G S for all f > 1, we have 


®Note that (291 remains valid if network coding operations are allowed. 


CJ G conv(5) for all T > 1. Since \U\ is finite, by the 
definition of lim inf, there exists a subsequence {u ■ CTfc}fe>i 
of the sequence {tt ■ Ct}t>1 such that 


min lim tt • Ct = min lim inf u ■ (32) 

uTzU k—^oo ^ uTzU T—^oo 


Since the set conv (5) C is closed and bounded, by the 
Heine-Borel theorem, it is compact. Hence any sequence in 
conv(5) has a converging sub-sequence. Thus, there exists a 
sub-sub-sequence {Cxk }*>i ^ conv(5) such that 

Ct^. ^ as t —> 00 . 

It follows that 


min tt • 

uGU 


if) 

(c) 


min lim tt • 

u£U i—^oo 

min lim tt • Cf 

uTz U k—^OD ^ 

min lim inf tt • Ct 

uGU T^oo 


min lim inf tt 

uGU T—>^oo 



where (a) uses the fact that if a;„ —^ a; then c ■ Xn ^ c ■ x 
for c, Xn, and x gM.\ I > 1; (b) follows that if the limit of 
a sequence {zn} exists then all subsequences {zn^} converge 
and limfc = lim„ (c) follows from Equation ( [3^ . This 
completes the proof of the lemma. ■ 

Combining Lemma with Eqn. ( [3T] i, we have that there 
exists a vector /3^ G conv(5) such that 


A* — e < min u ■ ff". 
ueu 


(33) 


Maximizing the right hand side of Eqn. 33 over all l3^ G 
conv(5), we have 


A* — e < max 

/3econv(5) 


( min u - (3] (34) 

y J 


Since the above inequality holds for any e > 0, by taking 
e \ 0 and expanding the dot product, we have 


A* 


< max 

j9Gconv(»S) 


min y Ce/3, 

a proper cut 

eGEu 


(35) 


B. Proof of Lemma 

Consider the wired network in Eig. |10(a)| where all edges 
have unit capacity and there is no interference constraint. Node 
a has total incoming capacity equal to two; thus, the broadcast 
capacity of the network is upper bounded by A* < 2. In fact, 
the network has two edge-disjoint spanning trees as shown in 
Eigures [l0(b)| and |l0(c)| We can achieve the broadcast capacity 
A* = 2 by routing odd-numbered and even-numbered packets 
along the trees 7i and T 2 , respectively. 

Consider a policy tt G Hin-order that ensures in-order delivery 
of packets to all nodes. Let Ri(t) be the number of distinct 
packets received by node i up to time t. Hence, node i receives 
the set of packets (1, 2,..., by time t due to in-order 

packet delivery. Consider the directed cycle a—>6—>c—>a in 
Eig. 10(a) The necessary condition for all links in the cycle to 
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(a) A wired network with a (b) Tree 7i 

directed cycle a b 
c a. 



Fig. 10; A wired network and its two edge-disjoint spanning 
trees that yield the broadcast capacity A* = 2. 


forward (non-duplicate) packets in slot t is Rait) > Rbit) > 
Rc{t) > Ra{t), which is infeasible. Thus, there must exist 
an idle link in the cycle at every slot. Define the indicator 
variable Xe{t) = 1 if link e is idle in slot t under policy tt, 
and Xe{t) — 0 otherwise. Since at least one link in the cycle 
is idle in every slot, we have 

(^) “f ^(fe,c) (^) “f ^(c,a)(0 — 

Taking a time average of the above inequality yields 

1 ^ 

y W + X(b,c){t) + X(^c,a}it)) > 1. 

t=l 

Taking a limsup at both sides, we obtain 


1 

limsup-^a;e(f) 

eG{(a,b).(t..c).(c,a)} ^ t=l 


> limsup E ^E Xe{t) > 1. 

eG{(o.b).(6.c),(c,o)} ^ i=l 

The above inequality implies that 

T 

1 , 1 

max limsup — > xJi) > -. (36) 

eG{(a.b),(b.c),(c,a)} T^oo T ^ 3 


Since the nodes {a, b, c} are symmetrically located (i.e., the 
graph obtained by permuting the nodes {a, b, c} is isomorphic 
to the original graph), without any loss of generality, we may 
assume that the link e = (a, b) attains the maximum in ( [3^ , 
i.e., 

T 

1 , 1 

lim sup - V X(a,b) {t)> o ■ (37) 

T—>-oo ^ 


Noting that Xe{t) = 1 if link e is idle in slot t and that node b 
receives packets from nodes r and a, we can upper bound 
Rb{T) by 


Rb{T) < ^ (1 

T 

<E(2 


It follows that 


X(a,b){t))- 


lim inf 

T —>-oo 


Rb{T) 

T 


< 2 


1 , 5 

limsup-Vx(a,b)(i) < X, 

T—)-oo ^ ._T ^ 


where the last inequality uses Thus, we have 

. . ,Rb{T) ^ 5 

min hm mt- < hm mt- < -, 

i&V T^ac T T-i-oo T 3 


which holds for all policies tt G Ilin.order- Taking the supremum 
over the policy class Ilin-order shows that the broadcast capacity 
Aj* jjrder subject to the in-order packet delivery constraint 
satisfies 

•^in-order = sup mm lim mf < - < 2 = A . 

7rGn.„.„,der T^OO T 3 


i.e., the network broadcast capacity is strictly reduced by in- 
order packet delivery in the cyclic network in Fig. |10(a)| 


C. Proof of Theorem 

We present the proof in three steps. First, using the dynam¬ 
ics of Xj{t) in d, we derive an expression of the one-slot 
Lyapunov drift using quadratic Lyapunov functions. Second, 
we design a useful stationary randomized policy that yields 
e-optimal broadcast throughput; this policy is useful to show 
that the system X (t), under the optimal broadcast policy tt*, is 
strongly stable for all arrival rates A < A*. Third, based on the 
above analysis, we show that the policy tt* is a throughput- 
optimal broadcast policy for any underlying network graph 
which is a DAG. 

Lemma 5. If we have 

Q{t + 1) < {Q{t) - + A{t) (38) 

where all the variables are non-negative and (x)^ = 

max{x, 0}, then 

Q\t + 1) - Q^{t) < + A2(f) + 2Q{t){A{t) - 

Proof: Squaring both sides of ( |38] l yields 

Q\t + l) 

< {{Qit) - + A\t) + 2A(t){Q{f) - gi{t))+ 

— {Q(f) ~ + 2A(t)(5(t), 

where we use the fact that x"^ > (a;+)^, Q(t) > 0, and /i(f) > 
0. Rearranging the above inequality finishes the proof. ■ 
Applying Lemma to the dynamics ( [l4| i of Xj (t) yields, 
for each node j r, 

X^it + 1)-X]{t) 
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<B{t)+2X,{t){Y^ 

(i') E (39) 

m^V keV 

where B{t) < + max{a‘^(t), < {a^{t) 

a{t) is the number of exogenous packet arrivals in a slot, and 
Mmax — maxeg E Ce is the maximum capacity of the links. We 
assume the arrival process a{t) has bounded second moments; 
thus, there exists a finite constant B > 0 such that E[i3(f)] < 
E(a2(f)) +2/i^^^ < B. 

We define the quadratic Lyapunov function L{X{t)) = 
Front ( [39] l, the one-slot Lyapunov drift A{X{t)) 

satisfies 


A{X{t)) A E[L{X{t + 1) - L{X{t)) I X(f)] 

= E[Y,{X^{t+l)-X]{t))\X{t)] 

j¥^r 

< B\V\ + 2^X,(f)E[ ^ ^ I X{t)] 

j^r mGV kGV 

= B\v\-2 

{iJ)^E kGKj{t) 

= B\V\-2 ^ E[fi,,{t)\X{t)]W,,{t), (40) 

where Kj (t) and Wij (t) are defined in ( fTS] ) and ( [T6] l, respec¬ 
tively. To emphasize that the evaluation of the inequality ( |40| ) 
depends on a control policy tt G If*, we rewrite ( |40l i as 

A-(X(f)) < B\V\ - 2 EK-W I X{t)] (41) 


Our optimal broadcast policy ir* is chosen to minimize the 
drift on the right-hand side of (|4T]i among all policies in If*. 

Next, we construct a randomized scheduling policy G 

If*. Let (3* G conv (5) be the vector that attains the outer 
bound on the broadcast capacity A* in Theorem [T] i.e.. 


(3* G arg max min 

/3Gconv(5) U: a proper cut 


eeEu 


From Caratheodory’s theorem ifTSll . there exist at most (|i?| -f 
1) link-activation vectors s; G 5 and the associated non¬ 
negative scalars {pi} with = 1’ ^rich that 

r = (42) 

1=1 

Hence, from Theorem [T] we have, 

A* < min CejS*. (43) 

U : a proper cut 

eGEu 


Suppose that the exogenous packet arrival rate A is strictly less 
than the broadcast capacity A*. There exists an e > 0 such that 
A -I- e < A*. From ( |43| ), we have 

A-I-e < min Ce/3*. (44) 

U: a proper cut ‘ ^ 

e^Eu 


For any network node v ^ r, consider the proper cuts Uy = 
V \ {u}. We have, from ( |44l l, that 


A + e < ^ Ce/3*, Vu ^ r. 
eeEu^ 


Since the underlying network topology Q = {V, E) is a DAG, 
there exists a topological ordering of the network nodes so 
that: (i) the nodes can be labelled serially as {ui,... ,n|y|}, 
where ui = r is the source node with no in-neighbours and 
v\v\ has no outgoing neighbours and {ii) all edges in E are 
directed from Vi ^ Vj, i < j ll20l : From ( |45] l, we define 
<71 G [0,1] for each node vi such that 

qi E Ce/3e = A-f e-^, ( = 2,...,|t/|. (46) 

Consider the randomized broadcast policy G If* work¬ 

ing as follows: (i) it selects the feasible link-activation vector 
s{t) = Si with probability pi in ( |42l ), I = 1,2,..., \E\ + 1, 
in every slot f; (ii) for each selected link e = {■,vi) of node 
vi such that Se{t) = 1, the link e is activated independently 
with probability qp, (iii) activated links are used to forward 
packets, subject to the constraints that define the policy class 
n* (i.e., in-order packet delivery and that a network node is 
only allowed to receive packets that have been received by 
all of its in-neighbors). Note that this randomized policy is 
independent of the state X{t). Since each network node j is 
relabeled as vi for some I, from ( |46l l we have, for each node 
j ^ r, the total expected incoming transmission rate satisfies 


5: E[7.r“wi^w]= E 

i:{i,j)&E i\{i,j)^E 

= qi Cel3l 


eeEu^^ 


(47) 


Equation ( |47l i shows that the randomized policy pro¬ 

vides each network node j ^ r with the total expected 
incoming capacity strictly larger than the packet arrival rate A 
via proper random link activations. According to the abuse of 
notation in at the source node r we have 

r rand . r X—> -^RAND . 

Ekr W I =E[ Y it)] = A- 

i:{i,r)GE i:{i,r)£E 

(48) 

From 63 and ( |48| ), if node i appears before node j in the 
aforementioned topological ordering, i.e., i = vi. < vi- — j 
for some k < Ij, then 


^ E[7tr“(t)i^(i)] 

k-.{k,i)^E 



k:{k,j)eE 


(49) 


The drift inequality ( |40l l holds for any policy tt G H*. Our 
broadcast policy tt* observes the system states X (t) and seek 
to minimize the drift at every slot. Comparing the actions taken 
by the policy tt* with those by the randomized policy 
in slot t in (|40ll, we have 


A-'{X{t))<B\V\-2 Y EK'{t)\X{t)]W.,{t) 

{i3)eE 

<i?|C|-2 ^ E[pf’'\t)\X{t)]W,,{t) 

{iJ)eE 


( 45 ) 
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= B\V\+2Y,X,it)( ^ I X(t)] 

j^r ^ m^V 

kev 

^50) 

j^r 

Note that = argmiiiiginQ) Qijit) for a given node j. Since 
node i’l is an in-neighbour of node j, must lie before j in 
any topological ordering of the DAG. Hence, the last inequality 
of ( |50l l follows directly from ( |49| l. Taking expectation in ( |50| ) 
with respect to X{t), we have 

E[LiX{t + 1))] - E[L(X(f))] < B\V\ - ^E||X(f)lli, 

where || • ||i is the £i-norm of a vector. Summing the above 
over t = 0,1, 2,... T — 1 yields 


T—1 

E[L(X(T))] -E[L(X(0))] < B\V\T-^ S2 

' ' t=o 

Dividing the above by 2Te/|y| and using L(X(t)) > 0, we 
have 

t^O ^ ^ 

Taking a lim sup of both sides yields 

limsup;^ y]] (51) 

T—^OO ■_p, 


which implies that all virtual-queues Xj(t) are strongly stable. 

Next, we show that the strong stability of the virtual 
queues Xj (t) implies that the policy tt* achieves the broadcast 
capacity A*, i.e., for all arrival rates A < A*, we have 


lim 

T—>-oo 


Rj{T) 

T 


= A, Vj. 


Equation d shows that the virtual queues Xj{t) have 
bounded departures (due to the finite link capacities). Thus, 
strong stability of Xj (t) implies that all virtual queues Xj (f) 
are rate stable Hi Theorem 2.8], i.e., Xj{T)/T = 

0, a.s. for all j. It follows that, 

lim — -= 0, w.p. 1 (52) 

T-J-OO T 


Now consider any node j ^ r in the network. We can construct 
a simple path a(r = Un ^ Un-i ■■■—?■ ui — j) from the 
source node r to the node j by running the following algorithm 
on the underlying graph Q{V,E). 

This algorithm chooses the parent of a node u in the path a 
as the one that has the least relative packet deficit as compared 
to u. Since the underlying graph G{V, E) is a connected DAG 
(i.e., there is a path from the source to every other node in 
the network), the above path construction algorithm always 
terminates with a path a{r j). The number of distinct 
packets received by node j up to time T can be written as 
a telescoping sum of relative packet deficits along the path a. 


Algorithm 3 r ^ j Path Construction Algorithm 

Require: Graph G(V^E), node j GV 
i ^ 1 
2: U, ^ j 

3: while Ui ^ r do 

4: Ui+i 3— argmin/jgin(„.) (5fciii(f); ties are broken arbi¬ 

trarily. 

5: i ^ i + 1 

6 : end while 


i.e., 

R,{T) = Ru,{T) 

n—1 

= Y, {Ru^ (T) - Ru.+1 (T)) + Ru^ (T) 

n—1 

= -Y^uAT)+Rr{T) 

n-1 T-1 

= -y]x„.(r)-f (53) 

i=l t=0 

where the third equality follows the observation that (see ( [T0| l) 

Xu.iT) = Qu,+^uAT) = Ru,+AT) - RuAT). 

Using Yl^iAi XuAA < < 1111 * and that Xj(f) are 

non-negative, we have, for each node j, 

T—1 T—1 

f E ^ E ^ ^ ^ E 

t—O j^r £—0 

Taking a limiting time average and the strong law of large 
numbers for the arrival process, we have 

lim = A, Vj. w.p. 1 

T—>00 1 

This concludes the proof 
D. Proof of Lemma 

We regard the DAG G as a wired network in which all links 
can be activated simultaneously. Theorem]^ and ([^ show that 
the broadcast capacity of the wired network G is 

A* = Adag = niin Y Ce = min Y Ce 

{/: a proper cut ' ^ 

eGEu 1 7- / eGEu^ 

= min din(n), (54) 

Dey\{r} 

where Uy = V \ {t;} is the proper cut that separates node v 
from the network, Ejj^ is the set of incoming links of node v, 
and the last equality follows that the maximizer in (|^ is the 
all-one vector f3 = 1 and that all links have unity capacity. 
Edmond’s Theorem Q states that the maximum number of 
disjoint spanning trees in the directed graph G is 

k* = min 'Y ^e- (55) 

U: a proper cut 

e^Eu 

Combining (|54ll and (|55]l completes the proof. 














